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Abstract 

In recent work, it was shown that velocity-dependent forces between paraUel funda- 
mental strings moving apart in a D- dimensional spacetime implied an expanding universe 
in D — 1-dimensional spacetime. Here we expand on this work to obtain exact solutions 
for various string/brane cosmological toy models. 
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1. Introduction 

In a recent paper [|l|, it was siiown tliat tlie velocity-dependent forces between parallel 
fundamental strings in D spacetime dimensions, with certain initial conditions, lead to 
an expanding universe in D — 1 dimensions. The findings were consistent with recent 
observations of an accelerating universe, and predict an asymptotially constant late 
time expansion rate. 

We start with the action S = Id + 82, where 

lD = i^j d^x^ge-^^ + A{dct^f - ^iff) (1.1) 
is the L>-dimensional string low-energy effective spacetime action and 

is the two-dimesional worldsheet sigma-model source action, qmn, Bmn and (j) are the 
spacetime sigma-model metric, antisymmetric tensor and dilaton, respectively, while 7^,^ 
is the worldsheet metric. H^, = dB2 and n is the string tension. The fundamental string 
solution to the combined action, representing stationary macroscopic strings parallel to 
the direction, is given by 

ds"^ = [-dt^ -f {dx^)"^) -f Sijdx'dx^, 

k (1-3) 

where n = D—4, = x'x^ and the indices i and j run through the L> — 2-dimensional space 
transverse to the string. The constant kn = 2K^Ti/nO„_|_i, where Ti = n is the tension of 
the string, equal to its mass/length, and O^+i is the volume of 5'"^"'"^, the n -|- 1-dimensional 
unit sphere. 

This solution can be extended to a multi-static string solution owing to the existence 
of a zero-force condition. This condition in turn arises from the cancellation between the 
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attractive gravitational and dilatational forces of exchange with the repulsive antisymmet- 
ric field exchange, and is based on the existence of supersymmetry and the saturation of a 
BPS bound |]. 

It was subsequently shown that, in addition to the zero static force, the leading order 
(0(f^)) velocity-dependent forces cancel for moving strings as well (see also 0). This 
result too is associated with the existence of higher supersymmetry . Following , it 
is straightforward to verify that the four-point amplitude corresponding to the scattering 
of two such fundamental string states approaches zero in the small velocity limit. This is 
identical to the result found for the a = -\/3 black holes, which also preserve half of the 
total spacetime symmetries, the maximum for such black hole, string or p-brane solutions. 

The Lagrangian for a test fundamental string moving in the background of a parallel 
source string is then given by []T|J^] 

C = -mh-^ i^^l - hx^ - l) , (1.4) 

where m is the mass of the string, = x'^Xi and the "■" represents a time derivative. 
It was shown in |]T| that the velocity-dependent force following from this Lagrangian is 
repulsive whenever the strings are moving away from each other, and this leads to a further 
separation of the strings. Since this type of interaction occurs for any two strings, if we 
start with any number of close, parallel strings initially moving apart in the transverse 
space, they will continue to do so indefinitely and will fuel an expanding universe in the 
D — 2-dimensional transverse space and therefore in the D — 1-dimensional spacetime 
orthogonal to the strings. For example, five- dimensional fundamental strings lead to an 
expanding universe in D = 4 spacetime dimensions. 

The toy model presented in consisted of a large number of fundamental strings 
initially very close to each other. Each pair of such strings interacts as above, so that 
an initial outward propagation of the strings tends to further push them apart in the 
transverse space. In a mean-field approximation, the effective force on each string was 
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approximated by that of a single, very large source fundamental string whose Noether 
charge k is equal to the total charge of all of the strings in the Z)-dimensional space. The 
distance r between the test string and the source string in this model then represents the 
approximate average position of the strings, and hence the size of the universe. The time 
dependence of r at both early and late times was determined for this expanding model. 

In this paper, we expand on the results of to obtain exact solutions for the radial 
position as a function of time for the mean-field approximation for D = 5 and D = 6 for 
the case of zero angular momentum. We also consider a spherically symmetric toy model 
and obtain similar results in the very early and late time limits as well, to give further 
support to the mean-field approximation results. 



2. Generalized p-Branes 

Before outlining these solutions, it is interesting to note that the Lagrangian ( |1.4| ) also 
arises whenever we consider the motion of a maximally supersymmetric p-brane moving in 
the background of a parallel, identical p-brane. For example, (|1.4|) is the same Lagrangian 
one obtains for a test fivebrane moving in the background of a parallel source fivebrane or 
for a DO-brane moving in the background of a source DO-brane (see, e.g., Q). This can 
be seen immediately either from supersymmetry, or through dimensional reduction 0. 
For the case of the fivebrane, for example, one replaces the two form Bqi with a six-form 
^012345 and proceeds in the same manner to obtain (|1.4[ ), where now m represents the 
mass of the test fivebrane moving in the background of a parallel source fivebrane. 

The relevant dimension is the number of transverse dimensions, given by n = D—p—3, 
since the harmonic function h = 1 + kn/r'^ depends only on n. Here 

kn = (2.1) 

where Tp, the tension of the p-brane, is equal to its mass/p-volume 0. Compactifying 
q < p dimensions, we can relate the Z)-dimensional Newton's constant Gd to the D — q- 
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dimensional Newton's constant G^-q via |^] 

GD = f^l = kI_^V, = Go-qV,, (2.2) 

where Vg is the compactified g- volume. Since Tp = m/Vp, it follows that 

= (2.3) 

which is just ( |2.1| ) for a p — g-brane. In particular, for g = p — 1, we recover the string 
formula. For q = p, we obtain the formula for DO-branes. As long as the longitudinal 
directions of the branes are held parallel, the dynamics are independent of p, the dimension 
of the branes, and depend only on the number of transverse directions. 

In what follows, we will consider parallel strings for simplicity, keeping in mind that 
we could equally well consider, say, DO-branes in one less dimension. 



3. Mean-Field Approximation 

For the case in which we replace the total repulsive force on a single string by an 
effective, large string at the origin, it was shown in [0 that 

where E is the constant total energy of the string, and where p = E/m is the ratio of the 
energy to the rest energy of the test string. We have set the angular momentum I = 0. 
Following, Chebyshev's Theorem ^, only the two cases of D = 5 (n = 1) and D = 6 {n = 2) 
can be integrated exactly. 



^ In order to integrate (|3.1| ) we need to evaluate integrals of the form J x"^(a+bx^)^dx (so called 
binomial differentials), where m, n,p are any rational numbers and a, 6-any constants. Chebyshev 
proved that integrals of this form can be expressed through algebraic, logarithmic and inverse 
circular functions in only three cases: 

1) p is an integer (positive, negative or zero) 

2) is an integer 

3) + p is an integer. 
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A straightforward integration for n = 1 yields 



p + 3\ / pr [r 7\ /y [r 7 / (p + 2)^ ft — t 



where a — -j^- For smaU r (or early time t), r c:::^ ^^/^i' "^^i^^ large r (or 

late time), r ~ a/ p{p + 2)/(p + ^ both in agreement with the findings of In this 
three-dimensional transverse space, r(t) represents the mean size of the universe in this toy 
model. Restoring factors of the speed of light c in (|2.1j ), it follows that ki = = 
where L is the length of each string and M is the mass of the source string, representing 
the effective total mass of the universe in the mean-field approximation. 

In it was claimed that r << fci and r >> fci corresponded to early and late times 
(relative to the current epoch) , respectively, in the expansion of the universe as implied in 
(p72|). Let us verify this assumption using estimates of cosmological parameters obtained 
in Ijl^. The current matter density p is very close to the critical density po = 3Hq/8tcG4, 
and the age of the universe to = Hq^ , where Hq is the Hubble's constant at present time. 
Ignoring numerical factors of 0(1), we take the current size of the universe ro ~ cto and 
its mass as M ~ pr^. It is then straightforward to show that ki is at most an order of 
magnitude less than tq. It follows that r « ki corresponds to much earlier times than 
present t « ki/c and r » ki corresponds to much later times t » ki/c. The model 
for n = 1 has the defect that the ratio of relative velocities to relative positions is not 
immediately a spatial constant, unless the spatial dimensions of the universe are restricted 
to two (see next section). 

For the more interesting case of n = 2, a straightforward integration yields 



^/p + 2\, r + Vr^^T~a- ^/a\ p(p + 2), , 

\p + lj \r + Vr^+a + y/al V {p + ly 



^ The constant energy E does not include the constant rest energy m. It is straightforward 
to show that, in the late time hmit, 7 = (l — Z?'^) ^^'^ = p + 1, so that the interaction energy is 
ultimately transformed into kinetic energy (= (7 — 1) m). 



where again a = and to is a constant. For small r, r c:^ rgexp while for large 
r we again find r ~ ^p{p + 2)/{p+l)H, both again in agreement with the findings of 
1^. In this four-dimensional transverse space, the three-dimensional universe may be 
regarded as an expanding spherical shell with radius r{t). A subtle point here arises as 
to the connection beween G5 and G4. In going from a five-dimensional universe to a 
four-dimensional one whose constant time slices consist of the expanding three-sphere, 
it follows that ~ G^rit). So G^ and G4 cannot both be constant. For constant 
G5, k2 = G^M (from (|2.1|) ) is also constant, but the four- dimensional Newton's 
law changes with time as the universe expands. The alternative picture is to demand a 
constant G4, but then allow for changing /c2, so that the mean-field model in this case 
should be thought of as an approximation to a cosmological p-brane solution with k2 a 
function of time. In either case, within the limits of these toy models, a straightforward 
calculation shows that /c2 ~ f^, where f is the current size of the universe. It again follows 
that the early (r << and late (r >> time limits are valid, as in the n = 1 case. 

Especially interesting features of the n = 2 case, other than allowing for a spatially constant 
Bubble's constant, are the infiationary expansion at early times and the asymptotically 
constant expansion rate for late times. This latter feature is generic to these string/brane 
models, since the velocity-dependent forces vanish at asymptotically large distances. 

4. Spherical Shell Model 

Now consider the following model of a string-seeded universe for both n = 1 and 
n = 2. N parallel, identical strings, with >> 1, are all located at the same distance 
R from the center of the transverse D — 2 dimensional space and move with the same, 
purely radial, velocity v = RR outward from the center. We would again like to determine 
R{t) and to compare our findings with the mean-field approximation. Before doing so, 
we note that such a model is consistent with the cosmological observation of a Bubble's 
constant. In the n = 1 (n = 2) case, the spatial universe consists of an expanding 2- 
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sphere (3-sphere). It is straightforward to show that the relative position is given by 
f2i = 1^21 1 = 1^2 ~ "H^ I = 2i?sin^/2, where 9 is the angle between the position vectors fi 
and r2. Similarly, the relative speed between two strings f 21 = |^2i| = 1^2 — "f^il = 2f sin 6*72, 
where vi and V2 are the velocities of the two strings. It follows that f2i/?'2i = f/-R is a 
constant over each sphere (or for a given time slice), representing the Hubble's constant 
for this model. 

The Hamiltonian for the system of test string moving in the source string background 



can be easily obtained from the Lagrangian (|1.4|) of this system 



H=^ ^ ' -1 , (4.1) 

For D = 5 [D = 6) this model is equivalent to a system of particles on the surface of a 
2-sphere (3-sphere) with two-particle energy of interaction. This interaction energy is just 
the difference between the conserved total energy of the 2-particle system given by ( [4.1J ) 
and the kinetic energy of the test string, since the source string is assumed to be stationary 
in the mean-field approximation. Note also that in the mean-field approximation (see also 
1^), the energy (|4.1| ) was taken to be a constant, which led to a solution for the motion of 



a single test string in the background of a much larger source string, which approximated 
the aggregate effect of the velocity-dependent forces of all the other strings. In the shell 
model, the interaction energies obtained from ( [4.1[ ) are not individually constant but must 
be added into a total energy for the system, which is then set equal to a constant total 
energy. The center of mass of the system remains at the center of the sphere. The 
interaction energy of one string in the background of another is then given by ^ 



^ It is convenient to use this velocity dependent form of the Hamiltonian. In the same way, 
one can easily obtain the conventional form H = H{r,p) by using the expression for the momenta 



Strictly speaking, we should use the more cumbersome relativistic form of the kinetic en- 
ergy. However, since most of the subsequent analysis involves the early time expansion, the 



where /ii2 = 1 + ^ and ri2 is the relative position of the strings. 

^12 

For D = 5 we assume for simphcity that the 1'^* string is located at the north pole of 
the 2-sphere with radius R. Then, by symmetry, the energy of interaction between the 1^* 
string and dN strings which are located inside the belt with azimuthal angles between 9 
and 9 + d9 is given by 



dEir,.fA 2 = dN 




1 ^ \ mr 
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1 - h^2rh I 2 



(4.3) 



where dN = f sin ^d9, ri2 = 2i2sin f , ri2 = 2i?sin f and hi2 = 1 + ^ 



In order to obtain the energy of interaction between 1*** string and all other strings we 
need to integrate ( [4.3| ) over 9 from 9 = to 9 = it. This can be done explicitly, but leads 
to a rather complicated, and not especially illuminating, expression. In order to make a 
connection with the early-time mean-field approximation, we first make the assumption 
that » 1. This means that the distance between any two strings ri2 << ki and our 
model can describe the system during the time when this condition holds, i.e. early times. 

Thus we replace hi2 — ^ in ( [4.3| ) . The integration over 9 is easy to perform, and 
we obtain for the energy of interaction between the 1*** string and all the other strings 

4mNR r 2 r . . 9 . ..all , 

= - — W (3a^ + 4a + 8) - 8] + - + - 1 (4.4) 



where 



a=— (4.5) 



Note that a is restricted to be in the domain < a < 1 for this model to be valid. 
By symmetry, the total energy of interaction of strings is 



non-relativistic approximation used for this model is valid. Furthermore, as we shall see later, the 
result for late times is not affected by making this simplification either. 



8 



Eint — NEiint 

and the total conserved energy of the system is 



(4.6) 



NEi=E = E,^t+N'^, (4.7) 

where Ei is the total energy of a single string. 

First assume a << 1 and expand the total interaction energy of the system ( [4.7|) in 
powers of a 

E^nt-^^N^a' + 0{a')]. (4.8) 
Thus up to 2^'^ order in a, Eq. (|4.7D takes the form 

3NR ^ R , 

where p = Ei/m = E/mN is again the ratio of the total energy (not including the rest 
energy) of each string to its rest energy. 



Solving the quadratic equation in ( [4.9|) for a, and using R « ki, it easily follows 



that the linear term in a in (|4.9|) may be dropped. It then follows that up to an 0{1) 



numerical factor, a ~ a / jt^ << 1. Since p is at least of 0(1), it follows that from a « 1 



that R » Dropping the 2"^'^ term in the r.h.s. of (0), solving (^) with respect to 
R and integrating we obtain 



,4.10) 

Another domain of interest is when a ~ 1. This corresponds to an even earlier time, 
since from ( fd.lOj ) it follows that a ~ so that an earlier time corresponds to larger a. 



For a — s> 1, the solution R = R{t) is easily obtained from the definition of a 
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(4.11) 

ki 



and is valid for R ~ ki/N, which can be seen from (|4.7|) with Eint given by ( [4.6| ) and 



Eiint by (|4.4|) with a ^ 1. Note that this quadratic expansion in time is consistent with 
the early-time approximation of the mean-field model. 

For large R {R >> ki), we can assume that h = 1 + ki/ R ^ 1, so that dH/dR = 0. 
From ([4.3| ), i\i.T\ ) it then follows that the constant energy Ei depends only on R and ^ 
so that R depends only on and p. So the radial velocity R is constant. Thus R oct for 
large R, again in agreement with the mean-field approximation. 

As mentioned above, for D = 6, the transverse motion of parallel strings is equivalent 
to the motion of particles with two particle energy of interaction given by ([4 .21) with 



hi2 = l + 4- (4-12) 



^2 



where ri2 = 2R sin ^ , ri2 = 2R sin ^ and x is the azimuthal angle, where we assume that 
the 1** string is located at the north pole of the 3-sphere. 

If we assume (as before) that the >> 1 strings are distributed homogeneously 
on the surface of the 3-sphere, then the number of strings located inside the belt with 
azimuthal angles between x ci-nd x + is 

dN = sin^xdx (4.13) 

TV 

and the energy of interaction between the 1** string (at the north pole) and dN strings 
inside the belt is given by the same expression (|4.3| ) with hi2 given by (|4.12| ) and dN by 



( [1.13|) . Assuming (as for D = 5) that >> 1 (i.e. R << ^fhi) and replacing h\2 
'"12 

we obtain 



^ This last statement is obviously also valid for the correct relativistic expression for the kinetic 
energy. 
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2Nm 



int 



sin xdx 



AR^ sin^ I 




1 I - 2i?2 sin^ ^ 



(4.14) 



integrating this expression over x from to tt we obtain the energy of interaction between 
the 1** string and aU other strings 



int 



2mNR'^ 




1 



mNR^ 



fl2 



(4.15) 



and conservation of energy condition can be written in the same form ( [4.7] ) as before with 
Ei^t given by ( |4.15| ) or 



1 



2NR'^ 



where 



k2R^ 
R^ 



(4.16) 



(4.17) 



where < 6 < 1. For R « ^Jh^N , it foUows from ( [4.16|) that 1 - 6 << 1. Alternatively, 
expanding ( 4.16|) in powers of 1 — 6 and keeping only the first nonvanishing term we easily 
obtain that condition 6^1 (or 1 — 6 << 1) leads to i? << \fk2fN . From the definition of 



h (|4.17|) we see that R ^ Rq exp in this case, again corresponding to the exponentially 
infiationary expansion in the D = 6 mean-field model. 

On the other hand, the condition << 1 or >> (again p is at least of 

0(1)) leads to 6 << 1, which can also be easily seen from ( [4.16 ). As before, expanding 
( [4.16| ) in powers of 6, keeping the lowest nonvanishing term and then integrating (in order 
to get R — R{t)) we obtain 



- J ^ , (4.18) 

For large R (i.e. R» v^) we can drop the 1 in Eq ( [4.12| ). In this case we obviously 
obtain the same result as for D = 5 (and generally any D): expansion with constant 

11 



radial speed, once more in agreement with the mean- field limit. We emphasize again that 
the nonrelativistic approximation for the kinetic energy in both cases does not affect the 
results for either early or late times. 

An interesting possibility in this case is that the moving strings in D = 6 lead to an 
expanding five-dimensional universe, in which an effective four-dimensional brane universe 
resides, following [^|. The asymptotic late time expansion rate is also intriguing, and may 
represent a testable prediction for this type of model. 

One possible advantage to the type of asymptotically fiat universe shown in these 



models is that, in contrast to a de Sitter universe, S-matrices would be well-defined 12 



At the same time, these models allow for an accelerating universe without assuming the 
existence of a cosmological constant. One can regard the changing acceleration as cor- 
responding to an effective cosmological "constant" which varies with time. For example, 
the -D = 6 (n = 2) model which has exponential growth in the early universe, has a con- 
stant effective cosmological constant, A ~ /c^^ O) which, however, is due entirely to 
the velocity-dependent forces between the strings/branes. At very late times, the effective 
cosmological constant is zero. It is then a straightforward but tedious exercise to determine 
the exact time dependence of the effective cosmological constant in these models. 

Further investigations of this type of model are clearly merited. More complicated and 
far more realistic models, possibly involving different species of branes could be considered. 
Furthermore, it would be interesting to go beyond the analytic, classical results obtained 
above, using a possible combination of numerical computations, quantum string effects and 
nonequilibrium thermodynamics. Nevertheless, it is likely that the leading order behaviour 
of the type of accelerating string/brane universe considered above is well-described in the 
classical approximation. Needless to say, these results await further verification and a 
better understanding of the underlying many-body interactions. 
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